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Abstract 


This  report  reviews  the  application  of  integral  equation  methods  that  are 
commonly  employed  to  study  interaction  of  EM  fields  with  two-dimensional  lossy 
anisotropic  objects.  The  study  partially  focuses  on  plane  wave  scattering  by 
arbitrarily  shaped  cylindrical  geometries  with  circular  and  rectangular  cross- 
sections.  In  addition  to  solid  material  objects,  investigation  of  the  bistatic  radar 
cross  section  of  several  metallic  circular  cylinders,  coated  with  a  single  layer  of 
anisotropic  material,  is  also  presented.  Finally,  radiation  effects  from  a  uniform  line 
source  through  a  lossy  anisotropic  shell  of  finite  thickness  are  also  investigated. 

The  report  includes  an  overview  of  the  theoretical  formulation  of  the  problem  in 
terms  of  the  combined-field  surface  integral  equations  (CFSIEs).  The  highligths  of 
the  Method  of  Moments  numerical  solution  to  CFSIEs  are  given,  and  numerical 
results  for  solid,  hollow,  and  coated  cylindrical  scatterers  are  presented.  The 
emphasis  in  the  discussion  of  the  numerical  results  is  placed  on  the  effects  of 
material  anisotropy  in  the  EM  response  of  the  scatterer. 


1 


Acknowledgment 


The  author  wishes  to  thank  James  C.  Logan,  Code  824,  of  NRaD  for  many 
helpful  discussions  regarding  the  numerical  results  appearing  in  this  report.  Thanks 
are  also  extended  to  Thinh  Q.  Ho  of  NRaD  for  providing  measured  values  of  complex 
dielectric  constants  of  several  composite  materials  used  in  this  study.  The  author  is 
also  grateful  to  Daniel  Tam,  Code  824,  of  NRaD  for  access  to  computer  resources 
that  were  used  to  generate  all  data  appearing  in  this  report. 


u 


Table  of  Contents 


Abstract . i 

Acknowledgment . ii 

1.0  Introduction .  i 

2.0  Integral  Equations  for  Anisotropic  Media . 3 

2.1  Formulation  of  Scattering  Problems .  3 

2.2  Formulation  of  Radiation  Problems .  7 

2.3  Method  of  Moments  (MOM)  Numerical  Solution .  9 

3.0  Numerical  Results .  12 

3.1  Plane  Wave  Scattering  From  Circular  Cylinders .  12 

3.2  Plane  Wave  Scattering  From  Cylindrical  Plates .  17 

3.3  Plane  Wave  Scattering  From  Coated  Circular  Cylinders .  19 

3.4  Radiation  From  a  Line  Source  Inside  a  Circular  Shell .  24 

4.0  Summary . 26 

5.0  References . 28 

Appendix:  Input/Output  File  Samples  of  MOM  Code  for 

Anisotropic  Cylindrical  Shells . A-i 


Section  1 


INTRODUCTION 


Numerical  modeling  of  EM  interaction  is  normally  performed  by  using  either 
differential  or  integral  equation  methods.  Both  techniques  can  be  implemented  to  solve 
problems  in  frequency  or  time  domain.  The  method  of  moments  (MOM)  approach 
(Harrington,  1968)  to  solving  integral  equations  has  matured  to  the  point  where  it  can 
be  used  to  solve  very  complex  problems.  In  the  past,  MOM  has  only  been  applied  to 
scattering  and  radiation  problems  involving  perfectly  conducting  or  isotropic  penetrable, 
lossy  or  losses  objects.  However,  many  materials  such  as,  for  example,  composites  that 
are  used  today  in  practical  applications  exhibit  anisotropy  properties. 

The  attractive  features  of  composite  media  include  the  combination  of  their 
mechanical  agility  and  correspondingly  desirable  electromagnetic  properties.  They  are 
much  lighter  than  the  metal  and  exhibit  comparable  strength.  Since  Navy  surface  ships 
are  built  primarily  from  the  metal,  they  are  beginning  to  reach  their  weight  limitations. 
In  order  to  consider  composite  materials  as  alternatives  to  the  metal  for  replacing 
portions  of  the  ship’s  superstructure,  detailed  studies  of  their  electromagnetic  properties 
are  needed.  Since  such  materials  are  anisotropic,  effective  numerical  tools  to  study  their 
interaction  with  electromagnetic  fields  are  required. 

To  account  for  these  effects,  several  integral  equation  formulations  for  scattering  and 
radiation  by  anisotropic  objects  have  been  developed  recently  (Graglia  et  al.,  1987; 
Monzon,  1988;  Beker  et  al.,  1989).  Integral  equations  are  usually  separated  into  those 
that  are  of  volumetric  or  surface  kind.  Although  the  volumetric  integral  equations 
(Graglia  et  al.,  1987)  are  very  general  and  can  be  used  to  analyze  scatterers  composed  of 
conductors  imbedded  in  anisotropic  or  inhomogenous  media,  they  require  solution  to 
matrix  equations  with  a  very  large  number  of  unknowns.  The  surface  integral  equations, 
on  the  other  hand,  are  best  suited  to  effectively  predict  the  scattering  response  from 
homogeneous  objects.  However,  partially  inhomogeneous  objects  can  be  analyzed 
employing  CFSIEs  as  well.  In  this  case,  the  entire  region  of  interest  must  be  subdivided 
into  multiple  locally  homogeneous  subregions  with  appropriate  boundary  conditions 
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enforced  to  connect  them  together.  Despite  the  added  complexity,  a  great  deal  of 
information  about  the  interaction  of  inhomogeneous  anisotropic  materials  with 
externally  generated  EM  fields  can  still  be  obtained  by  employing  them. 

This  report  is  aimed  at  examining  the  effects  that  material  anisotropy  has  on 
scattering  properties  of  arbitrarily  shaped  solid,  hollow,  and  coated  two-dimensional 
objects.  Numerical  results  are  provided  for  the  far-fields  of  lossless  as  well  as  lossy 
scatterers.  In  addition,  the  radiation  patterns  due  to  a  line  source  placed  inside  a  lossy 
anisotropic  circular  shell  which  are  functions  of  material  anisotropy  are  also  included. 
The  surface  integral  equations  for  all  aforementioned  geometries  are  solved  by  using  the 
standard  frequency-domain  two-dimensional  method  of  moments  codes  developed  at  the 
University  of  South  Carolina  (Beker  et  al.,  1989;  Beker,  March  1991;  and  Beker  et  al., 
1990). 

Although  most  practical  problems  are  three-dimensional,  a  great  deal  of  information 

on  the  electromagnetic  nature  of  composite  media  can  still  be  learned  from  studying  two- 

dimensional  geometries.  Radiation  of  a  line  source  placed  inside  a  shell  composed  of  lossy 

anisotropic  material  can  give  an  estimate  of  power  loss  due  to  reflections  and  dissipation 

* 

in  the  material.  At  the  same  time,  scattering  from  canonical  objects  such  plates  or  circles 
gives  an  idea  of  how  the  anisotropy  of  the  material  affects  the  RCS.  Consequently,  these 
studies  can  provide  sufficient  information  to  determine  preliminary  design  specifications 
for  the  composite  materials  that  can  be  potentially  used  for  parts  of  Navy  surface  ships. 
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Section  2 


INTEGRAL  EQUATIONS  FOR  ANISOTROPIC  MEDIA 


In  this  Section  the  highlights  of  two  kinds  of  surface  integral  equations  for  two- 
dimensional,  homogeneous  anisotropic  regions  are  presented.  The  formulation  to  several 
scattering  and  radiation  problems  is  given  and  its  method  of  moments  solution  is  briefly 
discussed.  Numerical  results  that  examine  effects  of  material  anisotropy  on  scattering 
and  radiation  from  arbitrarily  shaped  objects  are  provided  in  the  subsequent  Section. 

2.1  Formulation  of  Scattering  Problems 

Consider  a  homogeneous,  anisotropic,  and  arbitrarily  shaped  two-dimensional  region 
shown  in  Fig.  2-1.  The  medium  is  assumed  to  be  characterized  by  the  following  forms  of 
the  permittivity  and  permeability  tensors: 


y* 


o 


0 

0 

0 


(2-1) 


H oM  =  /*o 


Vy* 

0 


0 

0 


where  eQ  and  (l0  are  the  permittivity  and  permeability  of  free  space. 


(2-2) 


It  is  well  known  that  two-dimensional  fields  can  be  separated  into  their  transverse 
electric  ( TE )  and  magnetic  (TM)  parts  or  polarizations.  This  allows  for  the  duality  to  be 
applied  to  find  one  set  of  fields,  say  TE,  if  the  solution  for  the  other  set  (TM)  is  already 
known  (Beker  et  al.,  1989).  Therefore,  the  formulation  of  the  scattering  or  radiation 
problem  in  this  part  of  the  report  will  be  restricted  to  the  TM  polarization,  with  the 
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results  for  the  TE  part  obtained  via  duality.  In  this  report,  all  fields  and  sources  are 
assumed  to  have  a  time-harmonic  dependence  eJ<a  which  will  suppressed  throughout. 


Figure  2-1:  Segmented  contour  of  the  anisotropic  scatterer. 


According  to  (Beker  et  al.,  1990),  it  is  possible  to  constmct  the  integral  equations  for 
anisotropic  regions  in  terms  of  vector  potentials  only  or  to  use  a  mixed  (scalar  and 
vector)  potential  formulation.  Consequently,  the  combined  field  surface  integral 
equations  for  a  solid  (bulk)  scatterer  shown  in  Fig.  2-1,  exposed  to  the  TM  excitation,  can 
be  written  in  terms  of  vector  potentials  as 


Etl  —  jo)(Atl  +  At  2 )  —  2 


(vxe  MMvxS)' 
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e0£i 


(2-3) 
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where 
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kx  -  k0  elfdl 

(2-5) 

k]  =  kleay  =  (ji„nn  +fi%)  , 

(2-6) 

assuming  that  [iyi  =  -(1^. 

The  integral  equations  for  the  same  problem  can  also  be  recast  in  terms  of  mixed 
potentials,  which  can  be  achieved  with  the  help  of  the  Lorentz  Gauge  for  anisotropic 
media  (Beker  et  al.,  1990) 


(2-7) 


As  a  result,  equation  (2-3)  still  remains  the  same,  but  equation  (2-4)  now  becomes 


~h  xH;  =  -n  x  < 


jfl>(F1+FJ)  +  V(vrI+V'2)- 


VxA,  | 

VoVx  Vo 


(2-8) 


with  At,  F ,  and  y/  for  the  anisotropic  medium  given  by 


*  .=— M. 


(2-9) 
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(2-12) 


_  (x-xf  (y-yf 

"  V  H.  Hyy 

where  Jt,M,  and  pm  are  the  equivalent  electric  and  magnetic  currents  flowing  on  the 

surface  of  the  scatterer  and  magnetic  charge  density,  respectively.  The  magnetic  charge 
density  that  acts  as  a  source  for  pm  is  related  to  the  magnetic  current  density  M  by  the 

continuity  equation 

V-M  =  -jcopm  .  (2-13) 


Notice  that  equations  (2-8)  to  (2-12)  contain  expressions  for  the  potentials  in  the 
anisotropic  region.  These  same  equations  can  also  be  used  to  find  the  corresponding 
expressions  for  the  potentials  in  the  isotropic  region  by  setting  £a  =  £v  pa  =  p„  =  pu 
and  pv  =  jlyx  =  0.  At  this  point,  equations  (2-3)  and  (2-4)  or  (2-8)  represent  the  desired 
combined  field  surface  integral  equation  set  required  to  analyze  TM  scattering. 

» 

The  TM-polarized  incident  field  in  the  scattering  studies  is  assumed  to  be  a  plane 
wave  whose  electric  and  magnetic  field  components  are  given  by 


E‘  =  zE0e~U 

and 


(2-14) 


(2-15) 


where  tj  and  k  are  the  intrinsic  impedance  and  propagation  constant  of  Region  1, 
respectively.  The  corresponding  RCS  can  be  computed  once  the  equivalent  surface 
currents  have  been  determined.  The  z-component  of  the  scattered  E-field  is  calculated 
from 


E‘M  =  -|2|7i(p')W'«(W!)*'+i-|  (2-16) 

^  C  c 
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by  taking  the  large  argument  expansions  for  the  Hankel  functions.  In  the  above 
equation  ri  and  R  are  defined  by 


n  =  x  cos  <p'  +  y'  sin  <p' 


(2-17) 


(2-18) 


with  <p'  being  the  angle  that  the  unit  normal  makes  with  the  x-axis  at  the  point  of 
integration  on  the  contour,  and  where  p  corresponds  to  the  direction  of  far-field 
observation.  Finally,  the  RCS  can  be  computed  from  the  following  relation 


/ 

2\ 

2np 

E, 

K 

K 

y 

(2-19) 


with  Ez  obtained  from  equation  (2-16). 

» 

2.2  Formulation  of  Radiation  Problems 

To  illustrate  the  application  of  surface  integral  equations  to  radiation  problems,  the 
analysis  of  the  anisotropic  shell  (or  single-layer  radome)  illuminated  by  a  bne  source 
placed  inside  it  will  be  undertaken.  Figure  2-2  shows  the  geometry  under  consideration 
which  can  be  used  to  examine  properties  of  single-layer  two-dimensional  anisotropic 
radomes.  Although  such  structures  are  just  idealized  versions  of  the  actual  radome,  a 
significant  amount  of  information  about  the  influence  of  anisotropy  on  EM  fields 
radiating  through  a  material  layer  having  such  properties  can  still  be  obtained. 

In  general,  the  excitation  to  the  radome  could  be  delivered  in  the  form  of  an  external 
plane  wave  or  from  an  internally  placed  line  source.  For  the  plane  wave  incidence,  the 
formulation  of  the  scattering  problem  would  proceed  in  the  same  way  as  for  the  bulk 
object  considered  in  Section  2. 1.  However,  the  presence  of  an  additional  boundary  would 
require  for  the  combined  field  surface  integral  equations  (CFSIEs)  to  also  be  obtained  at 
the  newly  introduced  boundary  (Beker  et  al.,  1990).  Since  the  scattering  problem  for  the 
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radome  is  similar  to  the  radiation  problem,  at  least  in  its  integral  equation  formulation, 
the  surface  integral  equations  will  only  be  provided  for  the  latter  case. 


Figure  2-2:  A  two-dimensional  anisotropic  radome. 


The  boundary  conditions  (B.C.s)  on  the  tangential  field  components  at  boundary 
contours  Cj  and  C2  lead  to  the  following  set  of  CFSEEs  for  the  radome  shown  in  Fig.  2-2. 
At  the  inner  contour,  matching  E-  and  //-fields  yields: 


-£,;  —  jO)(AjX  +  Al2)  z 
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r 

jO)(Fl+F2)+V{y/i+ii/1)- 

^  /'oA  V*  Jj 

and  then  at  the  outer  contour  the  B.C.s  result  in 
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0  =  -nx{ 


jto(F3  +  F2)+V(yf  3+yr2)- 


VxA,  [^r(Vx^))| 
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(2-23) 


where  subscript  3  corresponds  to  Region  3,  which  for  this  problem  is  air.  Expressions  for 
the  vector  and  scalar  potentials  are  same  as  those  given  earlier  in  Section  2.1,  and 
incident  fields  are  those  of  the  line  source. 

For  this  problem,  the  TM-polarized  line  source  has  E  and  H  fields  that  are  given  by 
the  following  expressions  (Balanis,  1989) 


and 


(2-24) 


-  V  x  E' 
jo>M 


(2-25) 


where  It  and  f(p,<p)  are  the  line  source  current  and  taper  pattern  function,  respectively. 


2.3  Method  of  Moments  (MOM)  Numerical  Solution 

The  surface  integral  equations,  whether  they  are  cast  in  terms  of  vector  potentials 
only  or  they  are  stated  in  terms  of  mixed  potentials,  can  be  solved  numerically  by  the 
Method  of  Moments  (Harrington,  1968).  The  only  difference  between  these  two 
formulations  from  the  stand-point  of  the  numerical  implementation  is  the  complexity  of 
integrals  that  are  involved  in  filling  of  the  system  matrix.  In  the  Galerkin  procedure, 
both  expansion  and  testing  functions  are  taken  to  be  of  the  same  type,  and  for  all  results 
appearing  in  this  report  they  were  chosen  to  be  the  unit  pulse  functions. 

The  unknown  equivalent  electric  and  magnetic  surface  currents  ( Jt,M, )  are  first 
expanded  as  weighted  sums  of  pulse  basis  functions 
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(2-26) 
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W.(P')  =  2p“0',)M, 


(2-27) 


These  same  unit  pulse  functions  can  also  serve  as  the  testing  functions  for  E-  and  //-field 
equations 


'£(P)  =  I>*(PJ 


jn~l 


(2-28) 


2M 

<"(P)  =  ■  <2-29) 

m=W+l 

Since  the  scatterer  can  have  surface  discontinuities,  some  components  of  the  electric 
and  magnetic  fields  can  be  singular  at  these  points,  whether  the  objects  are  conducting 
(Jackson,  1975),  isotropic  (Bach  Andersen  et  al.,  1978;  Van  Bladel,  1985;  and  Marx, 
1990),  or  anisotropic  (Beker,  August  1991).  For  TM  scattering,  this  means  that  the 
tangential  //-field  (or  7X)  is  singular  at  the  surface  discontinuity.  On  the  other  hand,  the 

tangential  component  of  the  electric  field  (or  Mt)  is  continuous  at  that  same  location. 

To  properly  represent  such  field  behavior,  the  standard  procedure  is  to  stagger  the 
equivalent  surface  current  pulse  functions  as  shown  in  Fig.  2-1.  This  form  of  expansion 
builds  in  the  proper  field  behavior  numerically,  thus  avoiding  sampling  of  singular 
magnetic  fields  at  sharp  corners  or  numerically  violating  charge  conservation  of  the 
magnetic  charge  density. 

Once  expansion  and  testing  are  completed,  CFSIEs  are  converted  to  a  set  of  Unear 
matrix  equations,  which  for  the  bulk  scatterer  with  a  single  boundary,  take  on  the 
following  form: 
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(2-30) 


'[z~]  [^]}K]i=[[£„r 
>»]  [c]JL[w.]j  [[«i]  ■ 

The  individual  submatrix  elements  of  the  above  matrix  system  for  the  vector  potential 
formulation  (equations  2-3  and  2-4)  are  explicitly  defined  elsewhere  (Beker  et  al.,  1989) 
and,  therefore,  are  not  repeated  here. 

Although  for  anisotropic  scatterers  the  mixed  potential  formulation  (equations  2-3 
and  2-8)  has  been  implemented  numerically  (Beker  et  al,,  1990),  explicit  expressions  for 
the  matrix  elements  of  equation  (2-30)  do  not  appear  in  print.  However,  the  MOM 
approach  is  very  similar  to  that  used  in  CFSIEs  for  isotropic  scatterers  and  for  this 
reason  is  omitted  from  this  report  However,  an  interested  reader  is  referred  to 
(Umashankar,  et  al.,  1984),  among  other  references,  for  details  on  the  numerical 
procedure  required  to  reduce  the  CFSIEs  to  a  matrix  form  for  isotropic  scatterers. 


II 
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Section 

MOM  NUMERICAL  EXAMPLES 


This  Section  summarizes  numerical  results  obtained  for  anisotropic  scatterers.  These 
include  solid  circular  lossless  cylinders,  thin  lossy  plates,  coated  metallic  circular 
cylinders,  and  circular  shells.  The  combined  field  integral  equations,  presented  in  the 
previous  Section,  were  solved  using  the  Method  of  Moments.  The  MOM  codes,  developed 
at  the  University  of  South  Carolina,  were  written  in  FORTRAN.  All  data  were  generated 
on  an  IBM-compatible  486-50  MHz  PC,  using  version  4.01  of  Lahey  Fortran  Compiler. 
The  largest  matrix  size  of 480  x  480  was  required  for  anisotropic  shells  to  account  for  the 
unknowns  on  both  its  internal  and  external  contours.  The  solution  time  of  approximately 
10  minutes  was  needed  to  obtain  equivalent  surface  currents  and  bistatic  RCS  or 
radiation  pattern  for  a  single  frequency  point. 

3.1  Plane  Wave  Scattering  From  Circular  Cylinders 

* 

A  circular  anisotropic  cylinder,  illuminated  by  a  normally  incident  TE  polarized  plane 
wave,  is  shown  in  Fig.  3-1.  The  frequency  of  the  incident  field  is  assumed  to  be  300  MHz, 

and  the  electrical  size  of  the  object  is  k0a  =  Three  different  sets  of  anisotropic 

medium  parameters  were  chosen  to  illustrate  the  differences  in  scattering  properties  of 
such  cylinders  when  compared  to  their  isotropic  counterparts.  The  MOM  code  used  in 
computations  was  validated  in  earlier  studies  against  the  plane  wave  expansion  method 
(Fig.  3  of  Beker,  1989).  The  number  of  unknowns  for  all  circular  cylinder  studies  was  120 
for  each  equivalent  surface  current,  resulting  in  the  total  of  240  unknowns. 

As  mentioned  above,  three  different  sets  of  the  permittivity  tensor  elements  were 
considered:  (e„  =  2, e„  =  4),  (£„  =  4,en  =  2),  and  {ea  =  2,^  =  4,e^  =  -e^  =  2)  —  and  in 
all  cases  pLa  =  2.  Notice  from  the  integral  equations  presented  in  Section  2  that  the 
formulation  for  the  TE  polarization  involves  only  transverse  elements  of  [e]  and  only  the 
axial  element  of  [pj.  According  to  duality  (Beker  et  al.,  1990),  scattering  by  TM  polarized 
excitation  will  only  involve  transverse  elements  of  [p]  and  the  axial  element  of  [£}.  This 


12 


means  that  under  TM  excitation,  the  cylinder  will  look  isotropic  to  the  incident  field,  if 
the  material  is  nonmagnetic. 


Figure  3-1.  Geometry  of  a  circular  cylinder  and  a  thin  cylindrical  plate. 
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Figures  3-2  through  3-5  show  comparisons  between  the  scattering  properties  of 
isotropic  and  anisotropic  circular  cylinders.  In  this  case,  the  relative  permeability  of  the 
isotropic  cylinder  is  fir  =  2  and  the  axial  element  of  [ft]  for  the  anisotropic  scatterer  is 
[ln  =2.  The  intent  of  these  numerical  studies  is  to  illustrate  how  the  scattering 
properties  of  highly  symmetric  objects  such  as  circular  cylinders  can  be  affected  by 
anisotropy  of  the  material.  The  RCS  of  circular  anisotropic  cylinders  whose  permittivity 
tensors  are  diagonal,  namely  (£„  -2,eyy  =  4)  and  (£„  =  4,£^  -  2),  are  plotted  in  Fig.  3-2 
along  with  the  RCS  of  the  isotropic  scatterer  whose  permittivity  is  =  2.  As  can  be  seen, 
the  RCS  of  the  three  cylinders  differs  considerably,  especially  in  the  angular  range  of  90 
to  180  degrees.  This  graph  also  contains  a  very  interesting  feature.  More  specifically, 
when  the  scattering  angle  becomes  90  or  270  degrees,  the  RCS  of  all  three  cylinders  is 
same! 

In  the  next  example,  the  scattering  cross  section  of  the  isotropic  cylinder  with  £,.  =  4 

and  m  =  2  is  compared  to  the  same  anisotropic  cylinders  that  were  mentioned  above. 

Figure  3-3  displays  the  corresponding  RCS  curves,  indicating  the  degree  of  dissimilarity 

between  them.  Notice  that  in  this  case,  the  RCS  level  of  anisotropic  cylinder  with 

(£a  -2,en  =  4)  (case  (b))  is  somewhere  between  the  RCS  of  the  other  two  scatterers. 

* 

Final  examples  of  RCS  calculation  for  circular  cylinders  are  shown  in  Figs.  3-4  and  3- 
5.  They  emphasize  the  changes  in  scattering  characteristics  of  symmetric  objects  that 
are  associated  with  the  addition  of  the  off-diagonal  elements  to  the  permittivity  tensor. 
The  most  evident  on  the  RCS  as  =  -eyx  —  2  are  introduced,  is  a  dramatic  reduction  in 

the  field  level  scattered  from  the  cylinder.  In  fact,  the  reduction  is  by  more  than  10  dB, 
as  compared  to  both  isotropic  or  anisotropic  scatterers.  This  is  obviously  attributed  to  the 
presence  of  and  £yx.  The  net  effect  of  these  elements  is  to  reradiate  the  scattered  field 
in  different  directions,  rather  that  focusing  the  energy  in  one  direction  only. 

The  secondary  effect  on  the  RCS  due  to  introduction  of  £xy  and  £yi  is  the  destruction 

of  the  symmetry  of  the  scattered  field  about  the  direction  of  incidence.  This  indicates 
that  even  though  the  geometry  of  the  object  is  highly  symmetric,  the  shape  of  the 
scattered  field  pattern  can  be  controlled  independently  by  its  material  parameters. 
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Figure  3-2.  RCS  Comparison  for  isotropic  and  anisotropic  circular  cylin  ’ers.  Case  (a)  e,.  =  2,  pr 
=  2,  Case  (b)  £XI  =  2,£w  =  4 ;fis  =  2,  Case  (c)  e„  =  =  2;/ic  =  2. 


Figure  3-3.  RCS  Comparison  for  isotropic  and  anisotropic  circular  cylinders.  Case  (a)  ^  =  4,  Py. 
=  4,  Case  (b)  e„  =  2,en  =  4;/i0  =  2,  Case  (c)  £„  =  4,e„  =  2 =  2. 
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Figure  3-4.  RCS  Comparison  for  three  different  anisotropic  circular  cylinders.  Cases  (a),  (b),  and 
(c):  (£„  =  2,en  =  4,^  =  -eyx  =  2;/*fl  =  2),  (ea  =  2,ew  =  4;/^  =  2),  (£„  =  4,  ew  =  2;  na  =  2). 
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Figure  3-5.  RCS  Comparison  for  three  different  circular  cylinders.  Cases  (a),  (b),  and  (c): 
(£„  =  2,e„  =  4,c^  =  -eyz  =  2  ,fia  =  2),  (er  =  2,/ir  =2),  (er  =  4,n,  =  2). 


3.2  Plane  Wave  Scattering  From  Cylindrical  Plates 

To  examine  effects  of  material  anisotropy  on  objects  with  less  symmetric  geometries 
than  that  of  a  circular  cylinder,  the  scattering  proper  ss  of  a  thin  plate  were  also 
computed.  The  anisotropic  plate,  shown  in  Fig.  3-2,  is  lossless  and  is  characterized  by 
three  sets  of  material  parameters,  namely  (£„  —  2,  eyy  —  4),  (£„  =  4,e^  =  2),  and 

(£„  =  2,£w  =  4,6^  =  -e^  ~  2)  or  cases  (a),  (b),  and  (c),  respectively.  For  all  these  plates 

the  zz-element  of  permeability  tensor  was  assumed  to  be  equal  to  2.  The  excitation  to  the 
scatterer  is  a  TE-polarized  plane  wave  whose  frequency  is  300  MHz.  Two  directions  of 
incidence  were  considered:  one  broadside,  the  second  inclined  45  degrees  with  respect  to 
the  y-axis. 

The  radar  scattering  cross-sections  for  both  angles  of  incidence  arejiisplayed  in  Figs. 
3-6  and  3-7.  As  shown  in  Fig.  3-6,  for  broadside  incidence,  the  RCS  patterns  of  plates 
with  different  diagonal  permittivity  tensors  differ  significantly.  This  is  particularly 
evident  in  the  backs catte ring  direction  (<pt  =  180").  At  this  angle,  the  RCS  corresponding 
to  case  (b)  is  nearly  30  dB  lower  than  the  calculated  RCS  for  case  (a).  This  result  is 
attributed  to  the  incident  field,  which  at  tips  angle  has  a  y -component  only  and  senses  a 
higher  value  of  in  case  (a)  than  in  case  (b),  that  is,  £„  =  4  rather  than  e„  =  2. 

Note  that  for  the  plate,  the  addition  of  the  off-diagonal  elements  to  [e]  has  an 
opposite  effect  on  the  RCS  as  it  did  for  the  circular  cylinder.  Specifically,  introduction\pf 
=  -£yx  =  2  seems  to  raise  the  RCS  level  in  both  forward  and  backs cattering  directions 

(case  (c)  in  Fig.  3-6),  whereas  for  the  circular  cylinder  (Fig.  3-4)  the  opposite  was  true. 
One  possible  explanation  of  such  behavior  is  that  the  plate  is  electrically  thin  and  its 
scattering  characteristics  are  governed  primarily  by  £„.  In  fact,  as  can  be  seen  from  Fig. 

3-6,  curves  (a)  and  (c)  are  not  that  different  from  one  another,  except  the  asymmetry  and 
somewhat  higher  RCS  levels  for  case  (c). 

As  the  direction  of  incidence  is  skewed,  the  scattered  field  from  the  plate  changes 
considerably.  Under  such  conditions  the  incident  field,  which  now  has  bothx-  andy- 
components,  senses  every  element  of  [e]  equally.  Consequently,  as  illustrated  in  Fig.  3-7, 
there  are  no  drastic  differences  in  the  RCS  of  all  three  plates.  Observe  that  similarity 
between  cases  (b)  and  (c)  is  still  in  evidence,  despite  the  change  in  angle  of  incidence. 
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Figure  3-6.  RCS  Comparison  for  three  different  anisotropic  cylindrical  plates,  <p(  =  0°.  Case  (a) 
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Figure  3-7.  RCS  Comparison  for  three  different  anisotropic  cylindrical  plates,  <p.  =  —45°.  Case 
(a)  =  2,en  =4,  Case  (b)  e„  =  4,e„  =  2,  Case  (c)  £„  =  2,en  =  4,^  =  -en  =  2. 
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In  earlier  examples,  the  anisotropic  medium  was  assumed  to  be  characterized  by 
individual  tensor  elements  whose  numerical  values  were  quite  different  from  one 
another.  For  the  following  examples,  the  medium  parameters  correspond  to  those  of 
actual  composite  media  that  were  determined  from  the  measurements.  Numerical 
results  for  the  RCS  were  obtained  for  an  E-glass  plate  having  ea  and  £n  values  ranging 
from  5-  >0.1  to  5-  j 0.125  and  from  4. 5 ->0.09  to  4. 5 ->0.1125,  respectively.  Since  E- 
glass  is  non-magnetic,  both  fls  and  fJ.r  in  the  numerical  simulations  were  set  equal  to  1. 
The  excitation  to  the  plate  is  still  in  a  form  of  a  TE-polarized  planu  wave  with  the 
frequency  of  300  MHz,  striking  it  broadside. 

Figures  3-8  and  3-9  show  the  comparison  between  the  RCS  patterns  of  lossy  isotropic 
and  anisotropic  plates.  The  isotropic  scatterer,  in  this  case,  is  characterized  by  a  complex 
dielectric  constant  £r  =  5  -  >0.125.  It  is  apparent  that  the  small  degree  of  anisotropy  in 

E-glass  is  responsible  for  lowering  the  RCS  level  by  1  dB  over  the  entire  range  of  the 
scattering  angle.  Once  again,  the  scattering  properties  of  the  cylindrical  plate  seem  to  be 
governed  by  the  yy-element  of  the  permittivity  tensor.  Since  both  the  real  and  imaginary 
parts  of  £n  are  smaller  than  the  real  and  imaginary  parts  of  £, ,  both  forward  and 

backscattering  power  levels  of  the  anisotropic  plates  are  lower  than  those  of  their 

> 

isotropic  counterpart,  as  was  expected.  This  behavior  emphasizes  the  fact  that  even  low 
anisotropy  can  affect  the  response  of  the  scattering  object  to  the  incident  field,  especially 
to  its  polarization. 

3.3  Plane  Wave  Scattering  From  Coated  Circular  Cylinders 

The  data,  shown  in  Figs.  3-11  through  3-14,  concentrate  exclusively  on  the  effects  of 
anisotropy  on  the  RCS  of  conducting  cylinders  that  are  coated  with  a  single  layer  of 
lossless  anisotropic  material.  The  normalized  radius  of  the  conducting  core  is  kQb  =  1, 
and  the  radius  of  the  coating  layer  is  kQa  =  2  (see  Fig.  3-10).  As  in  all  the  previous  cases, 
the  coated  scatterers  are  also  excited  with  a  TE-polarized  plane  wave  whose  frequency  is 
300  MHz,  and  which  is  incident  at  <pt  =0°. 

Results  in  Figs.  3-11  and  3-12  convey  the  differences  in  RCS  plots  of  lossless 
anisotropic  (£a  =  5,  £n=  4.5)  and  isotropic  coatings  with  £r  being  4.5  or  5.  At  the  same 

time,  the  data  displayed  in  Figs.  3-13  and  3-14  show  the  comparison  of  RCS  patterns 
corresponding  to  the  anisotropic  coating  whose  diagonal  elements  of  [e]  are 
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Scattering  angle  in  degrees 

Figure  3-8.  RCS  Comparison  for  isotropic  and  anisotropic  cylindrical  plates,  <p,  =  0°.  Case  (a) 
e  =  5~ j  0.125,  Case  (b)  £a  =  5-j0.125,£w  =4.5-./0.09. 


Scattering  angle  in  degrees 

Figure  3-9.  RCS  Comparison  for  isotropic  and  anisotropic  cylindrical  plates,  <p;  =  0°.  Case  (a) 
£  =  5  -  ;0. 1 25 ,  Case  (b)  e„  =  5  -  j 0. 1,£„  =  4. 5 -  jO.  1125. 
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(£„  =  4.5,£y),  =  5)  and  those  of  the  isotropic  coating  layers  whose  £,  is  4.5  or  5  The 

scattered  field  levels,  for  the  most  part,  follow  the  expected  trends,  with  the  RCS  being 
higher  for  anisotropic  coatings  whose  en  is  greater  than  £,,  and  lower  when  is  less 
than  £, . 


Figure  3-10.  Geometry  of  a  coated  conducting  circular  cylinder. 

In  addition  to  the  expected  behavior,  numerical  results  for  coated  objects  have  other 
features.  For  example,  there  are  pattern  cross-overs  and  wide  angular  sectors  where  the 
RCS  of  isotropic  as  well  as  anisotropic  scatterers  are  identical.  The  greatest 
discrepancies  in  the  RCS  plots,  however,  are  confined  to  forward  and  backward 
scattering  directions,  as  well  as  to  the  null  locations.  Such  scattering  behavior  is  not  the 
same  as  was  seen  earlier  for  solid  isotropic  and  anisotropic  circular  cylinders  (Figs.  3-2  to 
3-5). 


Recall  that  for  the  solid  circular  scatterers,  there  was  no  correlation  between  the  RCS 
of  any  anisotropic  cylinders  with  the  RCS  of  any  isotropic  cylinders.  In  fact,  the  two  sets 
of  patterns  had  very  little  in  common.  The  same  can  no  longer  be  sai  1  for  loaded  circular 
cylinders.  This  is  most  evident  in  Figs.  3-13,  wherein  there  is  hardly  any  difference 
between  the  displayed  RCS  patterns.  This  is  a  direct  consequence  of  having  a  relatively 
large  conducting  core  that  both  anisotropically  and  isotropically  loaded  scatterers  have 
in  common. 


21 


Bistatic  RCS  in  dB 


3.4  Radiation  From  a  Line  Source  Inside  a  Circular  Shell 


The  effects  of  material  anisotropy  on  the  radiation  characteristics  of  a  uniform  line 
source  placed  inside  a  circular  shell  centered  at  the  coordinate  origin  are  now  presented. 
The  motivating  factor  behind  the  choice  of  a  circular  geometry  for  the  radiation  studies 
was  the  need  to  emphasize  effects  of  the  material  rather  than  the  shape  of  the  object  on 
its  radiation  properties.  The  shell  was  assumed  to  be  lossy,  while  its  internal  and 
external  radii  were  chosen  to  be  0.1  and  0.125  meters.  The  TE-polarized  source,  whose 
frequency  is  400  MHz,  is  located  at  the  center  and  is  assumed  to  be  radiating  uniformly 
in  all  directions. 

Numerical  results  displayed  in  Fig.  3-15  correspond  to  power  loss  estimates  as 
functions  of  the  anisotropy  in  the  conductivity  of  the  material.  The  two  curves  that 
appear  in  this  figure  were  calculated  by  varying  one  element  of  [a]  while  keeping  the 
other  one  fixed.  The  real  parts  of  ea  and  en  of  the  shell  are  5  and  4.5,  respectively.  An 

estimate  of  the  power  loss  due  to  dissipation  of  power  in  the  material  and  due  to 
reflection  back  into  the  core  of  the  structure  were  calculated  along  the  direction  of  the 
main  beam,  according  to  the  following  relation 

A P  =  i - ; - 1 - ! - - —  . 

|£r,([<x]=o)| 

Notice  that  the  power  loss  in  not  the  same  for  increasing  loss  tangent  values  in  x- 
and  y-direction.  The  loss  is  higher  for  changing  values  of  <7a  when  compared  to  same 

changes  in  a This  behavior  is  rather  surprising  since  the  electric  field  of  the  TE- 

exdtation  senses  both  complex  elements  £a  and  en  almost  equally. 

Finally,  the  corresponding  radiation  pattern  of  a  line  source  located  inside  the 
circular  radome  is  shown  in  Fig.  3-16,  which  contains  a  very  important  feature.  Despite 
the  highly  symmetric  geometry  and  uniformity  of  the  excitation,  the  calculated  radiation 
pattern  is  not  uniform !  This  is  a  direct  consequence  of  material  anisotropy,  even  though 
the  actual  difference  between  £a  and  en  is  small.  This  added  directivity  to  radiation 

pattern  can  be  increased  even  further  by  selecting  larger  differences  between  the  various 
tensor  elements. 
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Power  lost  along  the  main  beam  (%) 


Figure  3-16.  Radiation  pattern  of  a  uniform  line  source  placed  inside  anisotropic  shell. 
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Section  4 


SUMMARY 


In  this  report,  combined-field  surface  integral  equations  were  shown  to  be  a  viable 
tool  for  the  analysis  of  electromagnetic  field  interaction  with  anisotropic  materials.  The 
method  of  moments  solution  to  the  CFSIEs  was  obtained  and  used  to  generate 
numerical  results  to  illustrate  effects  of  material  anisotropy  of  scattering  and  radiation 
properties  of  arbitrarily  shaped  two-dimensional  objects.  Several  canonical  geometries 
such  as  solid  and  loaded  circular  cylinders  as  well  as  cylindrical  plates  were  considered 
which  can  be  used  effectively  to  investigate  EM  properties  of  composite  materials. 

Numerical  studies  indicate  that  material  anisotropy,  no  matter  how  small,  can  have 
a  significant  effect  on  the  electromagnetic  response  of  scatterers  composed  of  such  media. 
This  is  especially  true  of  composite  materials  which  are  known  to  be  anisotropic.  The 
computed  radar  cross-sections  of  objects  with  highly  symmetric  geometries  show  that 
depending  on  the  values  of  the  tensor  elements,  their  scattered  fields  can  be  changed 
considerably.  Even  asymmetry  in  the  far-field  pattern  can  be  introduced  with  addition  of 
the  off-diagonal  terms  to  the  permittivity  tensor.  The  asymmetry  seems  also  to  be 
associated  with  a  significant  reduction  in  the  RCS  resulting  from  redirection  of  the 
incident  field  by  the  constitutive  properties  of  the  material. 

Additional  studies  show  that  when  the  scatterer  is  electrically  thin  and  elongated  at 
the  same  time,  such  as  the  plate,  then  its  scattering  response  is  largely  controlled  by  a 
single  tensor  element  that  is  most  heavily  exposed  to  the  incident  field.  Under  these 
conditions,  the  anisotropic  scatterer  behaves  much  like  its  isotropic  counterpart  having 
scalar  permittivity  equal  to  that  of  the  dominant  tensor  element,  which  is  most 
prominently  sensed  by  the  excitation.  This  implies  that  electrically  thin  objects  are 
sensitive  to  the  polarization  of  the  incident  field. 

As  opposed  to  bulk  or  solid  objects,  conducting  objects  that  are  coated  with  a  thick 
layer  of  anisotropic  material  are  less  sensitive  to  its  tensor  characteristics  or  the 
polarization  of  the  excitation.  From  the  numerical  results,  it  appears  that  the  conducting 
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core  lessens  the  effects  of  anisotropy  by  not  allowing  for  the  field  to  totally  penetrate  the 
composite  structure. 

The  radiation  studies  reveal  additional  features  of  anisotropic  material  interaction 
with  electromagnetic  fields  that  are  otherwise  unobtainable  with  their  isotropic 
counterparts.  In  the  circular  shell  radome  example,  it  was  found  that  the  radiation 
characteristics  of  a  uniform  line  source,  placed  inside  this  radome,  are  affected  by  the 
anisotropic  material  properties  of  the  radome.  In  particular,  it  was  observed  that 
directivity  can  be  introduced  to  the  radiation  pattern  by  the  presense  of  the  anisotropy. 
This  is  a  very  important  feature  for  it  can  be  used  to  provide  additional  degrees  of 
freedom  to  the  design  of  shielding  structures  that  are  impossible  to  achieve  with 
conventional  isotropic  materials. 

In  addition  to  introducing  directivity  to  the  radiation  pattern,  anisotropy  can  also 
affect  the  power  loss  of  the  radome.  This  power  loss  is  defined  as  a  fraction  of  power 
transmitted  through  a  layer  of  lossy  and  lossless  material,  and  hence  is  a  sum  of 
transmitted  power  reduction  due  to  reflection  and  dissipation  loss.  If  the  conductivity  of 
the  material  is  also  characterized  by  a  tensor,  then  the  level  of  power  loss  becomes 
dependent  on  the  polarization  of  the  excitation.  It  was  found  that  a  combination  of 
factors  such  as  the  real  and  imaginary  parts  of  the  complex  permittivity  tensor  produce 
different  levels  of  power  loss  through  the  radome. 

From  the  results  presented  in  this  report,  it  can  be  safely  concluded  that  the  effects 
of  material  anisotropy  are  not  negligible  in  the  investigation  of  their  interaction  with 
electromagnetic  fields.  They  have  to  be  taken  into  account  in  any  numerical  methods 
used  to  simulate  EM  response  of  structures  with  such  material  properties.  Even  a  small 
degree  of  anisotropy  can  be  important,  especially  in  problems  involving  high  power 
electromagnetic  fields. 
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Appendix 


Input/Output  File  Samples  of  MOM  for  Anisotropic 
Cylindrical  Shells 


The  input  section  of  the  SHELL.FOR  file: 


PROGRAM  SHELL 


C . . . 

c 

C  This  program  computes  induced  surface  currents  and  radiated 
C  fields  of  a  dielectric  and/or  permeable  cylinder  coated  with  a 
C  single  LOSSY  layer  of  anisotropic  material.  The  composite 
C  scatterer  is  excited  by  a  TM-polarized  line  source  located 
C  inside  the  bulk  region. 

C 

C  Note:  The  magnitude  of  line  source  |Ezi|  =  1.0  V/m 
C  For  TE  polarization  use  duality:  Ez  ->  Hz, 

C  [eps]  ->  [mu] ,  and  (mul  ->  [eps]  . 

C . 

c 

C  Author:  Benjamin  Beker 

C  Dept,  of  Electrical  &  Computer  Engineering 

C  University  of  South  Carolina 

C  Columbia,  SC  29208 

C  (803)  777-3469 

C  E-mail:  beker@charlie.ece.scarolina.edu 

C . 


.  . . .  Main  Program  Input . . . . 


.  Relative  medium  constants  . 

. . . .  Index  1  corresponds  to  the  external  medium  and  index  2 
to  the  internal  dielectric  core.... 


FREQ 

=  4.0E+08 

EPS ( 1 ) 

=  1.0 

MU  ( 1 ) 

=  1.0 

EPS33 

=  CMPLX ( 

1.0  , 

0.0  ) 

MU11 

=  CMPLX ( 

5.0  , 

0.0  ) 

MU22 

*  CMPLX ( 

4.5  , 

-  0.005 

MU12 

=  CMPLX ( 

0.0  , 

0.0  ) 

EPS ( 2 ) 

=  1.0 

MU(2) 

=  1.0 

N(l)  - >  #  of  unknown  current  pulses  on  the  ext. 


A-l 


C  scatterer  for  each  of  the  currents 

c 

c  N (2 )  - >  #  of  unknown  current  pulses  on  the  int . 

C  scatterer  for  each  of  the  currents 

C 

C  NOTE:  neither  n(l)  or  n(2)  should  exceed  150 

C  unless  the  matrix  arrays  are  redimensioned 

C 

C  NCP(l)  - >  #  of  contour  points  on  the  ext.  scatterer 

C 

C  RA(1)  - >  radius  of  the  ext.  circle  in  meters 

C 

c  RA(2)  - >  radius  of  the  int.  circle  in  meters 

c 

C  NL  - >#of  contours  (currently  NL  =  2  which 

C  corresponds  to  a  shell) 

C 

c . 


N{  1 )  =  120 

N(2)  =  120 

RA ( 1 )  =  0.1254 
RA(2)  =  0.1 
NL  =2 


The  output  file  FOR018.DAT: 


TE-Scattering  Data 


Pulse 

Real  Part 

Imag .  Part 

Magnitude 

-  Mz  (Ext 

.  Contour)  - 

— 

1 

-0.0022386 

0.0030364 

1.4221663 

2 

-0.0022378 

0.0030355 

1.4216969 

3 

-0.0022361 

0.0030335 

1.4207293 

4 

-0.0022335 

0.0030307 

1.4193100 

5 

-0.0022302 

0.0030270 

1.4174374 

6 

-0.0022261 

0.0030224 

1.4151201 

7 

-0.0022213 

0.0030170 

1.4123945 

8 

-0.0022157 

0.0030108 

1.4092689 

9 

-0.0022096 

0.0030038 

1.4057910 

10 

-0.0022029 

0.0029962 

1.4019868 

55 

-0.0022261 

0.0030224 

1.4151261 

56 

-0.0022302 

0.0030270 

1.4174254 

57 

-0.0022335 

0.0030308 

1.4193202 

58 

-0.0022361 

0.0030335 

1.4207281 

59 

-0.0022378 

0.0030355 

1.4216964 

60 

-0.0022386 

0.0030364 

1.4221785 

61 

-0.0022386 

0.0030364 

1.4221652 

62 

-0.0022378 

0.0030355 

1.4217049 

63 

-0.0022360 

0.0030335 

1.4207244 

64 

-0.0022336 

0.0030308 

1.4193224 

65 

-0.0022302 

0.0030270 

1.4174256 

110 

-0.0021957 

0.0029880 

1.3978842 

A-2 


111 

112 

113 

114 

115 

116 

117 

118 

119 

120 

-0.0022029 

-0.0022096 

-0.0022157 

-0.0022212 

-0.0022261 

-0.0022302 

-0.0022335 

-0.0022360 

-0.0022378 

-0.0022386 

0.0029962 

0.0030038 

0.0030108 

0.0030170 

0.0030224 

0.0030270 

0.0030307 

0.0030336 

0.0030354 

0.0030364 

1.4019806 
1.4057915 
1.4092811 
1.4123819 
1.415125 4 
1.4174368 
1.4193102 
1.4207398 
1.4216809 
1.4221810 

121 

-  js  (Ext.  Contour)  - 

1.0492311  -0.6342636 

1.2260408 

122 

1.0490825 

-0.6343570 

1.2259620  • 

123 

1.0486403 

-0.6346370 

1.2257285 

124 

1.0479072 

-0.6350989 

1.2253408 

125 

1.0468915 

-0.6357400 

1.2248049 

126 

1.0456034 

-0.6365533 

1.2241269 

127 

1.0440555 

-0.6375290 

1.2233131 

128 

1.0422623 

-0.6386582 

1.2223728 

129 

1.0402439 

-0.6399286 

1.2213173 

130 

1.0380151 

-0.6413270 

1.2201540 

175 

1.0440561 

-0.6375304 

1.2233144 

176 

1.0456054 

. -0.6365552 

1.2241296 

177 

1.0468925 

-0.6357415 

1.2248065 

178 

1.0479070 

-0.6351005 

1.2253414 

179 

1.0486450 

-0.6346405 

1.2257345 

180 

1.0490789 

-0.6343566 

1.2259587 

181 

1.0492308 

-0.6342649 

1.2260413 

182 

1.0490873 

-0.634360^ 

1.2259680 

183 

1.0486403 

-0.6346380 

1.2257291 

184 

1.0479070 

-0.6351004 

1.2253414 

185 

1.0468932 

-0.6357417 

1.2248074 

230 

1.0330383 

-0.6444526 

1.2175744 

231 

1.0356061 

-0.6428420 

1.2189035 

232 

1.0380172 

-0.6413302 

1.2201574 

233 

1.0402443 

-0.6399315 

1.2213192 

234 

1.0422643 

-0.6386592 

1.2223749 

235 

1.0440576 

-0.6375302 

1.2233156 

236 

1.0456040 

-0.6365535 

1.2241275 

237 

1.0468916 

-0.6357406 

1.2248054 

238 

1.0479085 

-0.6351001 

1.2253426 

239 

1.0486411 

-0.6346372 

1.2257293 

240 

1.0490836 

-0.6343574 

1.2259631 

241 

-  Mz  (Int.  Contour)  - 

-0.0028280  0.0027539 

1.4881147 

242 

-0.0028276 

0.0027535 

1.4879004 

243 

-0.0028268 

0.0027526 

1.4874673 

244 

-0.0028257 

0.0027514 

1.4868546 

245 

-0.0028243 

0.0027498 

1.4860368 

246 

-0.0028225 

0.0027478 

1.4850240 

247 

-0.0028204 

0.0027454 

1.4838473 

248 

-0.0028181 

0.0027428 

1.4825174 

249 

-0.0028154 

0.0027399 

1.4810292 

250 

-0 . 0028126 

0.0027366 

1.4794073 

295 

-0 . 0028225 

0.0027478 

1.4850277 

296 

-0 . 0028243 

0.0027498 

1.4860350 

297 

-0.0028258 

0.0027514 

1 . 48b£593 

298 

-0.0028268 

0 . 0027526 

1.4874625 

299 

-0.0028276 

0.0027535 

1.4879123 

300 

-0 . 002828C 

0.0027539 

1.4881091 

301 

-0.0028280 

0.0027539 

1 .4881"20 

302 

-0.0028276 

0.0027535 

1.4879084 

303 

-0.0028269 

0.0027526 

1.4874744 

304 

-0.0028258 

0.0027514 

1.4868551 

305 

-0.0028243 

0.0027498 

1.4860313 

350 

-0.0028096 

0.0027331 

1.4776763 

351 

-0.0028126 

0.0027366 

1.4794093 

352 

-0.0028155 

0.0027398 

1.4810289 

353 

-0.0028181 

0.0027429 

1.4825237 

354 

-0.0028204 

0.0027455 

1.4838436 

355 

-0.0028225 

0.0027478 

1.4850290 

356 

-0.0028243 

0.0027498 

1.4860337 

357 

-0.0028257 

0.0027514 

1.4868461 

358 

-0.0028269 

0.0027526 

1.4874763 

359 

-0.0028276 

0.0027535 

1.4879031 

360 

-0.0028280 

0.0027539 

1.4881083 

-  Js  (Int. 

Contour)  - 

“  — 

361 

2.1365283 

0.2950120 

2.1567998 

362 

2.1365638 

0.2950518 

2.1568403 

363 

2.1366673 

0.2951703 

2.1569591 

364 

2.1368406 

0.2953671 

2.1571577 

365 

2.1370780 

0.295639^ 

2.1574302 

366 

2.1373813 

0.2959845 

2.1577780 

367 

2.1377437 

0.2963988 

2.1581938 

368 

2.1381588 

0.2968771 

2.1586707 

369 

2.1386247 

0.2974153 

2.1592062 

370 

2.1391389 

0.2980073 

2.1597972 

425 

2.1370809 

0.2956388 

2.1574330 

426 

2.1373799 

0.2959833 

2.1577764 

427 

2.1377499 

0.2963982 

2.1581998 

428 

2.1381590 

0.2968758 

2.1586707 

429 

2.1386256 

0.2974134 

2.1592069 

430 

2.1391463 

0.2980063 

2.1598043 

431 

2.1396866 

0.2986447 

2.1604276 

432 

2.1402729 

0.2993254 

2.1611023 

433 

2.1408827 

0.3000392 

2.1618054 

434 

2.1415083 

0.3007784 

2.1625276 

435 

2.1421559 

0.3015370 

2.1632745 

470 

2.1402714 

0.2993243 

2.1611009 

471 

2.1396878 

0.2986449 

2.1604288 

472 

2.1391366 

0.2980057 

2.1597946 

473 

2.1386251 

0.2974140 

2.1592064 

474 

2.1381588 

0.2968763 

2.1586707 

475 

2.1377399 

0.2963975 

2.1581898 

476 

2.1373787 

0.2959838 

2.1577752 

477 

2.1370776 

0.2956391 

2.1574297 

478 

2.1368382 

0.2953668 

2.1571553 

479 

2.1366663 

0.2951705 

2.1569581 

480 

2.1365628 

0.2950521 

2.1568394 

A*4 


Angle 

Scatt .  Field 

RCS 

RCS (dB) 

Norm.  Power 

360 . 00 

1.5325131 

1 . 1949408 

0.7734638 

3 . 9  669  4  63 

359.00 

1 . 5825242 

1 . 1949575 

0.7735247 

3 . 9670071 

358.00 

1.5825577 

1 . 1950082 

0 .7737056 

3.9871910 

357 . 00 

1.5826125 

1 . 1950909 

0.7740095 

3  .  9874921 

356 . 00 

1 . 5826898 

1 . 1952076 

0.7744335 

3.9879155 

355.00 

1.5827891 

1 . 1953576 

0.7749784 

3 .9884610 

175.00 

1.5827912 

1 . 1953608 

0.7749902 

3 .9884725 

174.00 

1.5829117 

1.1955429 

0.7756516 

3.9891341 

173.00 

1.5830539 

1.1957577 

0.7764319 

3.98993  42 

172 . 00 

1.5832180 

1.1960056 

0.7773319 

3.9908144 

371.00 

1.5834035 

1.1962858 

0.7783495 

3.9918318 

170.00 

1.5836093 

1 . 1965969 

0.7794787 

3.9929612 

169 .00 

1.5838367 

1.1969404 

0.7807255 

3.9942081 

168.00 

1.5840840 

1 . 1973144 

0.7820820 

3.9955645 

167.00 

1.5843515 

1 . 1977187 

0.7835486 

3  .9970310 

166.00 

1.5846393 

1 . 1981540 

0.7851261 

3.9986088 

165.00 

1.5849469 

1 . 1986190 

0.7868118 

4.0002942 

5.00 

1.5827889 

1.1953574 

0.7749777 

3.9884601 

4.00 

1.5826900 

1 . 1952080 

0.7744348 

3.9879172 

3 . 00 

1.5826129 

1.1950915 

0.7740115 

3.9874940 

2.00 

1.5825579 

1.1950085 

0.7737099 

3.9871922 

1.00 

1.5825242 

1 . 1949575 

0.7735247 

3.9870071 

0.00 

1.5825129 

1.1949404 

0.7734625 

3.9869452 
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